Rapidity scaling of multiplicity and flow in weakly and strongly interacting systems 
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We examine the "naturalness" of the scaling of multiplicity and elliptic flow V2 with rapidity in 
weakly and strongly interacting systems. We argue that multiplicity scaling is relatively straight- 
forward to incorporate in existing ansatze, and that this scaling is insensitive to the transport 
properties of the system. On the other hand, we argue that the observed scaling of elliptic flow is 
problematic to describe within a hydrodynamic model (the Knudsen number K <C 1), but arises 
more naturally within weakly interacting systems (where the Knudsen number ~ 1). We conclude 
by an overview of ways proposed to make weakly interacting systems compatible with the absolute 
value of elliptic flow, and by indicating experimental probes which could clarify these issues. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Nq 



I. INTRODUCTION 

The azimuthal anisotropy of mean particle momen- 
tum (parametrized by it's second Fourier component 1)2), 
thought of as originating from the azimuthal anisotropy 
in collective flow ( "elliptic flow" ) , has long been regarded 
as an important observable in heavy ion collisions. The 
main reasons for this is that elliptic flow has long been 
thought to be "self-quenching" [l|, 0]: The azimuthal 
pressure gradient extinguishes itself soon after the start 
of the hydrodynamic evolution, so the final V2 is insensi- 
tive to later stages of the fireball evolution and therefore 
allows us to probe the hottest, best thermalized, and pos- 
sibly deconfined phase. 

In addition, as has been shown in 3], elliptic flow is 
highly sensitive to viscosity. The presence of even a small 
but non-negligible viscosity, therefore, can in principle be 
detected by a careful analysis of v 2 data. 

Indeed, one of the most widely cited (in both the aca- 
demic and popular press) news coming out of the heavy 
ion community concerns the discovery, at the relativis- 
tic heavy ion collider ( RHIC ), of a "perfect fluid", 
also sometimes referred to as "sQGP" (strongly coupled 
Quark Gluon Plasma) d^. The evidence for this claim 
comes from the successful modeling of RHIC v 2 by boost- 
invariant hydrodynamics [Tol - [l5j . The scaling of v 2 ac- 
cording to the number of constituent quarks further sug- 
gests that the flow we are seeing is partonic in origin 

Going further in our understanding is hampered by 
the large number of "free" (or, to be more exact, poorly 
understood from first principles) parameters within the 
hydrodynamic model: While the equation of state is 
thought to be understood from lattice simulations, the 
behavior of shear and bulk viscosity is quantitatively not 
known around T c , where it is ex pec ted the temperature 
dependence could be non-trivial |2l| - |24| . The same goes 
for the large number of second-order transport coeffi- 
cients. While we have some understanding of the initial 
transverse density of the system (its dependence on the 
transverse coordinate is thought to be either a "Glauber" 



superposition of p-p collisions [25| or a partonic semi- 
classical "color glass" [26|, we do not as yet have 
control over the degree of transparency of the system, 
the amount of transverse flow created before thermaliza- 
tion (thought to be necessary to make the data agree 
with particle interferometry measurements [28|). or of 
the interplay between the "medium" and the surrounding 
"corona" of peripheral p-p collisions 29, 3fJ. A model in- 
corporating "all physics" , therefore, is expected to have 
a lot of correlated parameters which will be highly non- 
trivial to disentangle. 

A tool with the potential of overcoming these difficul- 
ties is scaling naturalness. Experiments have collected an 
extraordinary amount of flow data, encompassing a wide 
range of Energy (y/s) ,centrality (parametrized by num- 
ber of participants N part ), system size (mass number A 
of the nuclei), rapidity y and pseudorapidity 77, particle 
species and transverse momentum (j>t)- 

The experimental data collected allows us to "scan" 
observables dependence on variables relevant to the the- 
ory, and to see if the observable change when the same 
variable is obtained in different ways (for example, flow at 
mid-rapidity of a lower y/s collision compared with flow 
at the fragmentation region of a higher y/s collision, at 
the same multiplicity density dN/dy). 

This is an important test, because many theories have 
only one such relevant scaling variable (RSV) expected 
to drive the observable's value. Conversely, many the- 
ories have several RSVs, which combine to produce the 
observable's value. 

In the first case, one expects the observable to depend 
on only one scaling parameter in a straight-forward man- 
ner, no matter how the scaling parameter was arrived 
at (eg, flow should depend on dN/dy only, no matter 
weather dN/dy corresponds to mid-rapidity at low y/s 
or high rapidity at high y/s). In the second case, one 
expects the scaling to fail (flows for the same dN/dy but 
different rapidity /y/s will be different), because the RSVs 
required to produce the scanned parameter will generally 
be different as one scans in energy,system size, centrality 
etc. If this expectation is not fulfilled, complicating the 
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model (e.g., adding more transport coefficients to the hy- 
drodynamic model) is not likely to create a scaling where 
there was none, because the extra features generally also 
addRSVsB 

If small dimensionless quantities (SDQs) ctj can be con- 
structed from some of the RSVs, the above statement is 
mathematically equivalent to Taylor-expanding the flow 



observable around the SDQs [31 



(observable) ~ A + A-^ai + B-^a-i + O (aic^) + •■• (1) 

where Ai^o leading order, are independent of on. 

In case of hydrodynamics, the parameters relevant for 
the dynamics are expected to be the initial temperature 
and chemical potential (setting the speed of sound c s , the 
free path l m f p ),the initial system size, and the intrinsic 
scales of the theory. In practice, the chemical potential 
is thought to be irrelevant for higher energy collisions, 
and the intrinsic scales of QCD combine to produce one 
thermal scale T c , thought to coincide with a dip in the 
sound velocity and a change in the viscosity over entropy 
density ratio, rj/s (the dip and the change would be sharp 
for a phase transition, or smooth for a cross-over) . These 
relevant parameters can be accommodated in SDQs, but 

• more than one SDQ is necessary: At least Knudsen 
number K , eccentricity e, speed of sound c s 

• A dimensionful theoretical scale (T c ) appears in the 
relevant SDQs: c s and K depend on T,with a scale 
given by T c 

• some SDQs are purely intensive (the speed of 
sound) while others (the Knudsen number) mix in- 
tensive and extensive variables 

the last two points mean scanning is possible, given a 
wide enough space in \fs, centrality and rapidity. The 
first point means flow observables are not automati- 
cally expected to scale easily with the multiplicity, since 
several combinations of relevant variables (higher initial 
temperature vs larger system) can underlie the same mul- 
tiplicity, yet produce different dynamics. 



1 As an example of how this works, a very popular, and physically 
reasonable refinement to models describing heavy ion collisions 
is the "core-corona model", assuming that the bulk of the sys- 
tem, evolving collectively, is surrounded by a peripheral corona 
which is an incoherent superposition of p — p collisions [29l . kill . 
While one can implement this in many different ways, from the 
assumptions of the model it follows that more central collisions 
have more "core" and less "corona" , while collisions with a dif- 
ferent y/s have a similar percentage of "core" and "corona" but 
different dynamics for "core" and corona. Hence, if one com- 
pares data at different centralities and multiplicities but same 
dN/dy, one should not expect any kind of scaling of observables 
with dN/dy in models incorporating core and corona dynam- 
ics. Conversely, if such a scaling is found, it is an indication 
that "core" /" corona" separation does not play a big role in the 
physics 



These considerations can be promoted to a test of the 
hydrodynamic model, as several "simple" patterns have 
been found, centered around the flow at different ener- 
gies, system sizes and rapidities scaling with a few global 
dimensionless parameters (such as eccentricity and mul- 
tiplicity density). As we will show, in fact, such simple 
dependence is generally incompatible with an expansion 
such as Eq. [T] in the context of hydrodynamics with a 
QCD equation of state, but is more compatible with a 
weakly interacting system, where the Knudsen number is 
of order unity (each degree of freedom interacts only once 
in the system's lifetime). 

We shall illustrate these points with one observable not 
yet fully examined within viscous hydrodynamic models: 
the vi dependence on particle rapidity y = taxih~ p z /E 
or pseudorapidity r\ — tanh" 1 p z j \p\, where p z is the mo- 
mentum in the beam direction and, y — rj for p z ^> px, m, 
true for most particles away from mid-rapiditj0. 

The reason for the comparative lack of interest for vi 's 
dependence of rapidity is both numerical and physical: 
Numerically, introducing a rapidity dependence in the 
system makes the problem much harder to solve. Physi- 
cally, it was thou ght that at mid- rapidity Boost-invariant 
initial condition [32J is a good description. Since this is 
also the denser region, where the "hot" QGP phase will 
be most prominent, concentrating on this region and re- 
lying on Boost-invariance to shield it from fragmentation 
effects in the projectile and target rapidity is a priori a 
sensible approach. 

Experimental data, however, seems to show at best a 
very narrow central plateau for both the multiplicity dis- 
tribution [H, HI and elliptic flow |35|. While data 
for rapidity dependence close to mid-rapidity is some- 
what ambiguous (the PHOBOS collaboration reports a 
triangular shape |[, while the STAR and PHENIX col- 
laborations suggest a trapezium following the multiplic- 
ity distribution [y,0]), it is clear that t> 2 starts decreasing 
[j, Q far from what is naively thought of as the unther- 
malized fragmentation region. 

Even more tantalizing are the systematics in both en- 
ergy and rapidity of these observables: It appears that 
limiting fragmentation [36j holds as well for A-A colli- 
sions [37} as for more elementary systems. Multiplicity 
at mid-rapidity also seems to be "aware" of limiting frag- 
mentation, since it follows, with no "sudden" transitions, 
the trend EH 



dN 
dy 



N par t In yfs 



(2) 



where N part is the number of participants and \fs is the 
center of mass energy. 



2 Unfortunately, as in the common notation, r\ can mean both vis- 
cosity and pseudorapidity, s can mean both the square of the 
center of mass energy and the entropy density, while S is the nu- 
clear overlap area. The reader unfamiliar with the field should 
watch for the context in which these letters are used to under- 
stand their meaning 
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Even more puzzling, V2 for A-A systems follows these 
trends as well, both w.r.t. to scaling with pseudo-rapidity 
[H| and at mid-rapidity [3lj . 

Fig. [T] summarizes these findings. It is the purpose of 
the current work to assess these experimental observa- 
tions in the light of current models used to understand 
RHIC data. We should underline that we are analyzing 
the scaling rather than the absolute values of these ob- 
servables. Detailed studies of the latter have been made 
at top RHIC energies [38l-l42|. and it was found that 
models combining ideal hydrodynamics and hadronic cas- 
cades reproduce multiplicity pretty well, and do a reason- 
able through not perfect job of reproducing vi- Scaling 
of these observables, within a theory containing intensive 
dimensionful parameters, has not however been quanti- 
tatively tested, and could hold the key to see whether 
data-model disagreements seen in |38l-l42j are compatible 
with experimental systematic errors, can be fixed with 
minor model refinements (coronae, rescattering, a more 
realistic initial state,...) or whether they indicate a more 
profound misunderstanding of the system under consid- 
eration. 

Section [TT| aims to demonstrate that scaling of dN/ dy 
(both in the limiting fragmentation region) and at mid- 
rapidity, is generally compatible with expectations from 
high-energy QCD, and that it is rather insensitive to 
the phase of the system at thermalization. The reader 
should not treat this section as a demonstration of a rig- 
orous model, but as an argument that a rigorous model 
in agreement with what we know from high energy QCD 
can be developed, and as an outline of the requirements 
such a model needs to satisfy to describe the scaling. 

The main section of this work, Section Mil demon- 
strates that the situation is rather different for the lim- 
iting fragmentation of elliptic flow, since, unlike dN/dy, 
elliptic flow should be sensitive to the phase of the sys- 
tem at thermalization, and that therefore the observed 
fragmentation of vi is not natural within hydrodynamics 

SectionlTvl argues that, in fact, such scaling arises more 
naturally in a weakly interacting theory, and section [V] 
outlines experimental probes capable of clarifying the sit- 
uation and closes with concluding remarks. 

II. SCALING OF dN/dy 

A. Limiting fragmentation 

While an exhaustive explanation for limiting fragmen- 
tation [36| still evades us, it is a plausible feature of many 
models easily justified within the QCD paradigm, in par- 
ticular in the models based around the "wounded nu- 
cleon" or "wounded quark" concept 43] (essentially, a 
superposition of independent collisions each producing 
elongated objects in rapidity). Without loss of gener- 
ality, we shall concentrate on the Brodsky-Gunion-Kuhn 
(BGK) [3,[45[ model in its more generic form (Fig[3](a)), 
to illustrate why limiting fragmentation is plausible, and 



to suggest why the logarithmic dependence of the mul- 
tiplicity density at mid-rapidity can be plausibly under- 
stood as being related to limiting fragmentation. While 
the BGK model is based on a generic scenario in terms 
of parton-parton collisions, its ansatz is easily visualized 
within the flux tube ("QCD string") picture. It should 
be noted that the Color Glass coordinate [2(| HtJ ini- 
tial state follows this ansatz, with the crucial addition of 
transverse dynamics dictated by saturation. The ansatz 
in [44| is however more general, and therefore likely to be 
applicable to regions in y/s and centrality beyond those 
where the assumptions behind [2(| are reasonable. 

In this picture, each collision between a projectile and 
target nucleon produces energetic partons, distributed 
uniformly in rapidity throughout the forward-traveling 
projectile and the backward-traveling target. These par- 
tons act as "string ends" ,on average for 
p{ymin) ~ Npart strings between the parton and the 
target, and p(y m ax) ~ Npart strings between the par- 

TP 

ton and the projectile (N't are the number of partic- 
ipants from, respectively, the target and the projectile). 
Medium partons are produced through fragmentation of 
these strings, uniformly in rapidity within each string. 
Partons from fragmenting strings are initially flowing 
with Boost-invariant flow (v z — z/t so the spacetime ra- 
pidity 

tanh- 1 = y 

tdyn 

(where Zdymtdyn are the formation spacetime points) , 
as expected in part on g eneration from longitudinally ex- 
panding flux tubes |46l | and confirmed by measurements 
such as net baryon number in rapidity [35] (showing that 
high xf — Pz/E valence partons tend to pass through 
the collision region with very little change in momen- 
tum). However, parton density p is not invariant with y, 
but rather follows from the sum of independent fragmen- 
tations of these flux tubes. 

It is easy to see 0, EE] that the resulting multiplic- 
ity distribution is a linear extrapolation from target (at 
Umax) to projectile (at y min = ~y max in the center of 
mass frame) 

Pbgk(v)~ A -By (3) 

where 

D j\j projectile jsjtarget 

_ = 2e Hp (4) 

A ^y-max 

and B depends weakly on y/s (g^j < and scales as 
the number of participants N part . 

All we did here is illustrate the features a theory should 
possess to describe limiting fragmentation in a natural 
way. While a detailed quantitative explanation of how 
these features arise in QCD evades us to this day, it is 
generally thought they are not unreasonable: As Fig. [5] 
(a) shows, it is natural to expect B to simply count Flux 
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FIG. 1. (color online) The experimentally observed multiplicity rapidity distribution (a), the multiplicity at midrapidity as a 
function of y/s (b) , and the rapidity dependence of V2 (c) . (a) is usually known as limiting fragmentation and (c) is referred to 
as limiting fragmentation of elliptic flow. All experimental plots from [341 ] . with the dashed line indicating the scaling trend 



Tubes (be they the old-fashioned electric flux tubes [46[ 
or the magnetic CGC version [13]), and these flux tubes 
to stretch between the liberated parton and the projectile 
or target rapidity. It is therefore not surprising that B is 
only weakly dependent on y/s 

In the projectile=target limit the BGK picture tends to 
the boost-invariant limit which, even at RHIC, does not 
seem to arise in experimental data [34l 35]. Furthermore, 
in this picture, it is not natural to reproduce the scaling 
Pbgk(u — 0) [= A] ~ lny^s seen experimentally. The 
next subsection discusses possible ways to fix this. 



B. scaling with y/s 

The most straight-forward interpretation of the 
4?7 ~ Np ar t In y/s scaling is that the limiting frag- 

mentation region reaches until mid-rapidity (the ini- 
tial partonic distribution is a triangle, and any central 
plateau is dominated by thermal smearing rather than 
longitudinal flow) . Since yu m ~ In y/s, assuming that 
the slopes of dN/dy distributions are constant with re- 
spect to y/s, the logarithmic scaling of the mid-rapidity 
multiplicity density would automatically follow. 

This straight-forward interpretation, however, is not 
trivial to reconcile with long-held beliefs regarding high- 
energy QCD used to develop ansatze such as BGK. Lim- 
iting fragmentation is only natural in a system with lit- 
tle stopping power, since it implies that, in the frame 
at rest with respect to the target, the projectile be- 
comes a Lorentz-invariant pancake while the target is 
extended [36[ . That is not the case in the center of mass 
frame, where both projectile and target should be equally 
Lorentz-contracted to a size much smaller than their rest- 
frame size. 

Furthermore, limiting fragmentation follows naturally 
if distributions created by each parton-parton collision 
are invariant in rapidity, as in the traditional BGK pic- 
ture (the linear rapidity dependence in non-central col- 



lisions comes from combinatorics of the strings created 
from the target and the projectile, not from microscopic 
dynamics). 

This, in turn, is thought to be based on asymptotic 
freedom, which suppresses multi-GeV momentum trans- 
fers necessary to shift degrees of freedom around the re- 
quired rapidity range. RHIC experiments have shattered 
our understanding of asymptotic freedom applied to 
multi-particle processes, and stimulated the development 
of several theoretical scenarios (strong semi-classical 
fields |48|, semi-classical small- a; physics [26|, strongly 
coupled quantum fields tractable by Gauge/string meth- 
ods [4!| and so on) where liberated partons interact 
"strongly". This raises the possibility that the stopping 
of color charges deviates from the naive BGK-like value. 

It is not immediately clear, however, how large stop- 
ping at mid-rapidity could produce limiting fragmenta- 
tion away from mid-rapidity. One needs not just to 
break boost-invariance, but to extend the set-up lead- 
ing to Eq. I .SI II so that the rapidity narrow trapezium is 
centered around mid-rapidity, and not around the pro- 
jectile/target local transverse parton density. 

For instance, Landau hydrodynamics (50| appears to 
exhibit approximate limiting fragmentation close to the 
edges of the rapidity distribution, although corrections, 
even neglecting transverse flow and using the confor- 
mal ideal gas EoS rather than the QCD one, are larger 
than the deviations seen in the data [51]. But the lim- 
iting fragmentation would not extend to anywhere near 
mid-rapidity, since Landau hydrodynamics at y <C yu m 
becomes indistinguishable from Bjorken hydrodynamics 
[H after, typically, a few T^ ia} ~ O (1) /y/s 0. Con- 
sequently, dN/dy at mid-rapidity in Landau hydrody- 
namics is not logarithmic wrt y/s (Hoj unless the initial 
temperature is unnaturally adjusted to produce the log- 
arithmic dependence. 

Similarly, weakly coupled transport models such as 
can produce rapidity distributions without plateaus, but 
generally break limiting fragmentation. This is not sur- 
prising, since limiting fragmentation is not natural in a 
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FIG. 2. (color online) The "usual" BGK picture discussed in section Hi A\ (sl) and the modified BGK picture discussed in section 
IIIBK b). In the usual picture, the partons of the projectile and target (denoted by red and blue filled circles) act as origin points 
of color flux tubes (denoted as green dashed lines). Partons from each flux tube are generated in a boost invariant configuration, 
but, after summing over all flux tubes, parton density is sensitive to the ratio of wounded nucleons in the projectile and the 
target, according to a linear extrapolation. The phenomenological modification in the (b) can be understood in terms of a third 
dominant source of color flux tubes, located at the spacetime locus of the collision between the projectile and target (denoted 
by the larger hollow magenta circle) 



model with multiple partonic scatterings. 

The simplest way to account for the data is to try to 
tinker with the BGK ansatz so that the "upright trapez- 
iums" seen in Fig. [5] (a) become triangles described 
in the first paragraph of this section, shown in Fig. [2] 
(b). One sketch of a way in which this could be done 
within the existing QCD language is to allow for enough 
(~ Apart + Nj>art) color charge stopping at the center of 
mass frame to act as a source of color strings along 
with the ends. Partons created in the initial scatter- 
ing (each scattering distributes them uniformly across 
rapidity, just like in the BGK picture) can then be joined 
by strings to the projectile, target, and the central colli- 
sion region. If the density of at the projectile and target 



3 These slow color charges are sometimes called, and might well be, 
"wee partons" l2(j| . For the discussion here, however, we shall 
just refer to them as "color charges" or "string ends", as they 
might not be partons at all, but some kind of collective excitation 
in strongly coupled dynamics a la [49| . All we need is that they 
carry color 



"string ends" 
and in the middle it is 



min,max J 



N, 



T,P 
part 



P 



Urain Umax \ \ AjP 

— j ^ iv p ar i r J*p a 



Np ar t + N part ^ P {Umin, Umax) 



then, following the ansatz used to derive Eqs. 13141 of each 
liberated parton producing a string connecting either the 
projectile, the target, or (most likely) the overlap region. 
If strings fragment to a boost-invariant distribution, we 
get 



p(y) ~ A PiT - B p ,t \y-y c 



where 



lie 



^projectile j\jtarget 

-^projectile ^ -^target 



(5) 



(0) 



in this ansatz, the initial rapidity distribution becomes 
not a boost-invariant distribution but a "triangle" (Fig. 
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[3J (a)) , which thermal smearing can reduce to the ob- 
served near-Gaussian (Fig. [3] (b)). Furthermore, since 
Umax ~ ln-y/s, and the slopes from the center of mass 
to the projectile/target are, respectively, N^ t (and in- 
dependent of y/s). It is easy to see that, for the linear 
extrapolation to be valid, 



A - (N^ art + N^ art ) InVJ 



Npart In ' 



(7) 



Essentially, the modification of the BGK picture we have 
made, together with the requirement of limiting fragmen- 
tation forces the scaling seen in Eq. [2] at mid-rapidity. 
Note that in the highly asymmetric limit (jp — A, or A — A 
at large transverse radius) this picture tends to the BGK 
one, since the center of mass also moves with large ra- 
pidity . 

Fig. El (a) summarize the initial parton distribu- 
tion at different energies expected from the modified 
BGK model, (b) simulates thermal smearing within the 
Bjorken co-moving frame, through a plot of f(y — y^, y^) 
defined as 



f(y,VL) 



VL 



VL 



dy' ( 1 - ) cxp [- cosh (y 1 



VL) 



(8) 

for two representative values of yi,. As can be seen, lim- 
iting fragmentation and the mid-rapidity abundance fol- 
low the parton distributions very well even after thermal 
smearing, while the overall rapidity shape looks very sim- 
ilar to the experimental one: The smeared triangle looks 
like a Gaussian. A modification of the triangle into a 
narrow trapezium, suggested by the plateau reported in 
[H, |j| in both dN/ dy and V2 will not modify this conclu- 
sion, provided the rapidity width of the tip is not much 
larger than one unit. 

The initial distribution of string ends in the above sce- 
nario is "three-pronged" , with peaks in th distribution at 
Xf ~ ±1 (mostly valence quarks, as can be seen from net 
baryon measurements [35j |) and a third peak at xp = 
(with width <C 1 . In the standard BGK picture this third 
peak would not be there). Strings would then stretch be- 
tween either of these peaks, and the newly created par- 
tons (evenly distributed in rapidity). 

It is not at all clear, and not the aim of this work, 
to specify which dynamics gives such a distribution of 
color charges. The question to answer here is: Why 
should string ends be created from strongly interacting 
(not necessarily non-perturbative!) dynamics at the over- 
lap point in configuration space, but only from usual 
parton-parton scattering and fragmentation elsewhere? 
We only note that, in spacetime rapidity, the loci of the 
system where color charges initially reside are exactly at 
the ones corresponding to these peaks: The lightcones 
of the target and projectile, and the spacetime region 
where the target and projectile overlap, at the center of 
mass rapidity. Thus, the distribution of color charges we 
conjecture is not unnatural if strong short-range (cut-off 
at a scale < Aqq D ) many-body dynamics at the overlap 
point is present. The gradients of these for high energy 



collisions are very large f~ y/sNp^j , so, provided the 

strong many body dynamics is cut off at large distances, 
color charges between \xp\ <C 1 and \xf\ ~ 1 (the par- 
tons between the three peaks, which function as other 
ends of the string) will only be produced by pQCD scat- 
terings, which, up to LHC energies, are subleading to soft 
particle emission. Strings completely disconnected from 
all three peaks would therefore be less likely. 

Once again, the preceding reasoning should be re- 
garded as a sketch of a way in which both dN/dy ~ 
Npart In \fs dependence and limiting fragmentation arise 
within the framework of QCD (partons, strings, etc.), or 
of a way to initialize initial conditions for a dynami- 
cal model (hydrodynamics, transport) which respects the 
scaling seen in experiment. It is not a model, since we do 
not specify which form of initial dynamics best satisfies 
Eq. I5I6I71 and because the microscopic dynamics inter- 
polating between initial partons and later strings is left 
undetermined. 

But we note that a model which does not has lit- 
tle chance of describing both limiting fragmentation and 
dN/dy ~ In y/s, as required to agree with experimen- 
tal data when a wide range of \/s (O (10 — 100) GeV) is 
scanned. 



C. dN/dy and the phases of QCD 

It is important to note that here we discussed the ini- 
tial formation state of the partons (the start of dynamics 
Tdyn), not the initial local equilibrium state of the system 
(the start of hydrodynamics r eq , subsequent to string for- 
mation, decay of strings into partons, and parton equili- 
bration). This is why we treated the system as partonic 
throughout the rapidity region. 

The properties of the medium at r eq need to be de- 
termined through the initial temperature, obtained the 
QCD equation of state to link the density to temperature 
and chemical potential. Thus, there will be a region in 
rapidity (\y\ < y c ) where, at r = r eg , the bulk of the 
initially equilibrated p will be carried by partons, and 
another region in rapidity where it will be carried by 
hadrons (put in in a different way, a region where the 
initial state will be a quark-gluon plasma, and a region 
where it is a hadron gas. For some energies the first re- 
gion dominates, for others the whole system has little or 
no QGP component). 

To roughly estimate the rapidity at which, on average 
over the fireball "slice" , the switch happens, we have used 
the ideal 3-flavor QCD equation of state to calculate the 
initial temperature assuming boost- invariant flow (not 
full boost-invariance) [32[ 



T — 

eq — 



dN 



4.9 



1/3 



dy Teq N*l 3 rt fm* 



■p(- 



eq) 



Nl/3 



(9) 



where g = 20.8 is the appropriate degeneracy constant 
for QCD. The transition point is around where T eq ~ 




FIG. 3. (color online) (a): The parton density in the modified BGK picture, (b): The final hadron distribution function, 
including thermal smearing. 
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FIG. 4. (color online) (a): The initial temperature dependence on rapidity for different y/s, assuming the multiplicity dis- 
tribution shown in Fig[T](a) (experimental data from [341] ). the Bjorken initial state given by Eq. [9] and r eq ~ T~ l . The 
symbol shapes and colors follow the coding of Fig. Q](b). Panels (b) and (c) show, respectively, the same calculation assuming 
T eq ~ an< i r eq ~ 1/ hi -^/s . The black ellipse in (a) indicates the point where effects due to the phase transition are 

expected to be most important, while red ellipses in (b),(c) indicate where limiting fragmentation would be expected to break 
down 



T c ,where T c ~ 180 MeV is the critical temperature for 
the phase transition [53| . Fig. H] (a) shows the result of 
this exercise: The critical rapidity for top RHIC energies 
is y — 3.5 — 4, corresponding to a critical y/s (where that 
region is at mid-rapidity) of \fs ~ 3 GeV (or 18 GeV 
fixed target, intriguingly close to the "kink, horn" and 
"step" observations |54j ) . It should of course be remem- 
bered that this is a very rough approximation, but not an 
unreasonable one considering that the ansatze described 
in the previous subsection mean that the pre-cquilibrium 
medium (at ~ Td yn ) is partonic throughout the rapidity 
range. 

The estimate shown in Fig. @] (a) was made assum- 
ing r eq ~ T~ ~ 1/^, constant w.r.t. yfs. This is not 
what is usually assumed: By the uncertainty principle, 
T eq should scale as 1 / y/s, or at least 1 / In y/s if saturation 
effects dominate. It should however be noted that any 
change of T eq in energy and system size would generally 
require an unnatural fine-tuning to reproduce the limit- 
ing fragmentation, as panels (b),(c) of the figure show. 
As Fig. 2] shows, the "transition" region is well within 
experimental acceptance, and hence any physical effect 



of it should be evident from a rapidity scan. 



To understand what this effect might be, we note that 
both nearly free streaming evolution and close-to-ideal 
liquid expansion conserve entropy. Therefore, the final 
dN/dy is relatively insensitive to whether the system was 
deconfined or not at equilibration. In other words, if at 
a given y and yfs there is a critical rapidity y c where 
a "perfect liquid" quark-gluon plasma changes into a 
"lousy liquid" hadron gas, we would not see "anything 
special" when measuring the dN/dy at y c , as entropy 
conservation would ensure final hadronic multiplicities 
would mirror initial partonic ones both in a long-lived 
sQGP (K <C 1) and a short-lived collection of nearly 
free-steaming hadrons (K ~ 1). Flow observables, and 
in particular vi , should however be very sensitive to such 
changes. The next section will demonstrate these points 
explicitly. 
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ratio rj/s with temperature (N c denotes the number of colors and A the 't Hooft coupling constant). Red ellipses indicate the 
temperature regime where one could expect scaling violations to occur. 



III. NATURALNESS OF SCALING OF 
ELLIPTIC FLOW IN HYDRODYNAMICS 

Assuming the rapidity dependence of flow observables 
is "encoded" in the initial density (and associated inten- 
sive properties: T,rj/s etc.) rather than in the trans- 
verse size (S ~ A 2 / 3 at all y), and assuming subse- 
quent evolution is local in y (a natural assumption in the 
weakly interacting limit, less so in the strongly interact- 
ing limit, since rarefaction waves can travel across rapid- 
ity [HI, 56] ) , we should expect the scaling of V2 to follow 
the scaling of dN/dy provided the intensive properties 
are either invariant with p(T eq ), or change monotonically 
with p(r eq ), throughout the rapidity range. However, in a 
cross-over from a QGP to a hadron gas, intensive proper- 
ties of the system should not change monotonically with 
p{T eq )- _ 

In Fig. [5l we summarize the expected changes: The 
viscosity to entropy ratio rj/s is expected to jump from 
the relatively high value of the hadron gas {rj/s ~ AT? 
in a gas of mesons and glueballs, where N c is the num- 
ber of colors [13]) to the low value of strongly interact- 
ing QGP (rj/s ~ O(0.1)iV c [H), and then slowly in- 
crease to the asymptotically free weakly coupled QGP 
(n/s ~ £2_{^ _2 ) ^ci where A is the 't Hooft coupling 
constant [59]). In addition, the speed of sound is expected 
to have a dip in the cross-over re gion , whose depth is at 
the moment not well determined [60, \QM ■ 

In a wide variety of models, elliptic flow V2 depends 
on the eccentricity e, and should be approximately pro- 
portional to it. This essentially follows from Taylor- 
expanding the solution of whatever dynamical equation 
V2 obeys in e, since e is small and dimensionless, and since 
by symmetry with no e there is no elliptic flow. 

We also know that V2 /e decreases if viscosity is turned 
on, i.e. if K increases. Hence, it is quite natural that, as 
suggested in [33, EH| 
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FIG. 6. (color online) The V2 dependence on rapidity 
given initial conditions reproducing limiting fragmentation, 
and subsequent hydrodynamic evolution. The superimposed 
dashed line shows the pre-equilibrium (r = Tdyn) partonic 
density from Fig. [3] (a) 



where K is the Knudsen number and Kq ~ O (1) is a pa- 
rameter specific to the theory. Furthermore, the trans- 
verse Knudsen number at a given mean free path l m fp 
is 



dN 



ImfpS dy 



(11) 



This formula assumes just boost-invariant flow, as well as 
a a time-scale [62| r„2 = \/~S/c s for the building up of V2, 
where (~ N^ t , Not to be confused with the center 
of mass energy y/s or the entropy density s) is the initial 
transverse size of the system (which, as we argued earlier, 
is independent of rapidity). It should be noted that going 
beyond this rough approximation for t V 2 worsens scaling, 
since [HH3] T v2 rises and saturates with increasing density 
due to the self- quenching of elliptic flow. 

The derivation of Eqs. [TU] and [TT] follows straight- 
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forwardly [62j from density formula 32 1 

1 dN 

P' 



Sr eq dy 



(12) 



(where S is the initial transverse surface and r eq ~ l m f P 
) and Taylor expanding 

We further remember that — ., , depends on the 

e \ideal 

equation of state, i.e. on the speed of sound. By a lead- 
ing order expansion argument, and remembering that the 
asymptotic expansion speed of a Godunov-type hydrody- 
namic shock wave ~ c s [E3, EH HH > it can be seen that 
— -j 7 ~ c s (numerical simulations lend credence to this 
scaling, see [62j). We also remember that l m f P ~ iS-. 

Putting everything together, and neglecting the differ- 
ence between y and the pseudo-rapidity rj (small in the 
fragmentation region away from mid-rapidity), we get 
that 



^ ^ c s (r e q ) 



part 



c s r\ 



Ts 



(13) 



we believe that when T > T c 77/s -C 1, c s ~ 1/V3, when 
T ~ T c c s -C 1/V3 and 77/s is at a minimum, and when 
T < T c c s goes back to a value not too different from 
l/v3 but 77/s increases to > 1. (Fig. [5]). 

Additionally, close to mid-rapidity the plasma lifetime 
should 3> r„2, so v 2 saturates and becomes independent 
of initial density (Eq. Q2] over-predicts 7j 2 /e). In the 
less dense region, however, the plasma lifetime < r v %, so 
v 2 should be approximately proportional to the initial 
density (Eq. [13] is a good approximation). 

On the other hand, N part should be independent of ra- 
pidity and pseudorapidity while the initial T(r eg ) should 
smoothly change as given by Eq. |H1 

We immediately see (Fig. [B]) that the scaling seen in 
Fig. Q] is n °t compatible with a modified BGK initial 
condition, or indeed any initial condition without an un- 
physically finely tuned correlation between the size of the 
system and intensive parameters 31]. In the supposedly 
long-lived ideal fluid mid-rapidity region, v 2 /e should be 
considerably flatter than dN/ dy due to the self-quenching 
of v 2 0- At the critical rapidity where T eq ~ T c , v 2 /e 
should dip due to the dip in the speed of sound, and in 
the fragmentation regions where To < T eq v 2 /e should 
go down more rapidly than dN / dy due to the rise in 
r)/s. The rapidity at which v 2 vanishes should in general 
be different from the rapidity at which dN/dy does, due 
to v 2 additional dependence on 77/ s and system lifetime. 
These considerations are qualitative, and based on some 
simplifying assumptions (c s and 77/s do not change over 
the timescale t v2 , and stay constant significantly above 



and below T c ). Because it is a scaling argument, however, 
more realistic dependence of c s and 77/s on temperature 
naturally leads to more scaling violation. 

While no detailed hydrodynamical studies of energy 
scans with 3D hydrodynamics (ideal or viscous) have, 
so far, been performed, the studies of 7J2 with pseudo- 
rapidity at the top RHIC energy [38ri42T ] confirm the in- 
sights of the present work. 

These simulations, while describing v 2 at mid-rapidity 
acceptably, have qualitative features absent from the 
data, such as fast changes in gradient of v 2 w.r.t. pseu- 
dorapidity (Fig. 10 of [39j], Fig 3 of [40]). Some choices 
of initial conditions minimize these differences ([38[ Fig 
3 and 4) but they never seem to disappear, as is natural 
due to the existence of intrinsic scales set by the phase 
transition and freeze-out conditions. 

These features are expected to become stronger if 
a temperature dependent viscosity (a la Fig. [S] (b)) 
is added to these models, and especially when non- 
equilibrium parameters, such as pre-thcrmalization flow 
(necessary for particle interferometry data) are included. 
Experimental data, however, has no trace of these 
changes in gradient. 

For one energy, this statement would only be as good 
as the non-negligible error bars of the experimental mea- 
surement. The limiting fragmentation of v%, however, 
would be spoiled had these variations of gradients been 
there. Given that this limiting fragmentation is seen, 
spanning y^i = 20 — 200 GeV [34J , the evidence that v 2 
drops uniformly with rapidity, with a gradient indepen- 
dent of *Js, seems to be solid. 

An obvious way to make this scaling more natural is 
to assume that, rather than 77/s being constant, we have 



V 

— ~ K 
S \ S 



T; 5 



(14) 



where k does not change across T c . Such behavior, how- 
ever, contradicts the intuition from kinetic theory (where 
77/s - Tl mfp - O (A~ 2 ) - O (T )). Furthermore, this 
solution is itself unnatural, since the equation of state ex- 
periences a cross-over and degrees of freedom change. For 
example, for SU(N C ) rj/s ~ N® above T c and ~ N^ below 
T c . Finally, this solution requires that the "partonic tri- 
angle" initial condition is also appropriate at lower (AGS) 
energies, where the initial temperature is according to the 
Bjorken formula < T c even at mid-rapidity. 

The observation of v 2 limiting fragmentation, there- 
fore, poses a formidable challenge to the hydrodynamic 
model. 



IV. SCALING OF ELLIPTIC FLOW IN THE 
WEAKLY INTERACTING LIMIT 



4 Note that, as mentioned in the introduction, the existence and 
size of the mid-rapidity plateau of V2 is currently controversial. 
Even if it exists, however, it can not be much larger than one 
unit of rapidity 



The discussion in the previous section motivates us to 
try to see if such a scaling could arise in weakly interacting 
systems, where collective effects play a smaller role and 
dynamics is determined just by the dilution and, rates of 
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scattering and mean fields. The advantage of this is that 
the latter two could, in certain circumstances, be less 
sensitive to whether the microscopic degrees of freedom 
at equilibration are partons or hadrons. 

In fact, a scaling of vi with density can be naturally 
derived by solving the Boltzmann equation with an ellip- 
soidal profile in initial transverse density, as was done in 
(65j . This gives rise naturally to a scaling relation of the 
form 

v 2 (av) dN 

where S is once again the transverse overlap region, and 
(av) ~ lmfp/p is the interaction cross-section between 
the degrees of freedom times their average relative veloc- 
ity. This scaling is surprisingly similar to that derived 
in Eq. QTj except that av is replaced by (a) c s ■ In a 
dilute system there is no reason, of course, to suppose v 
is related to c s rather than to some initial intrinsic mo- 
mentum. For example, for massless particles v = 1 while 
c s = 1/V3. 

Together with the modified BGK scaling discussed pre- 
viously, such an ansatz will naturally lead to a scaling of 
the type observed, provided that (av) does not vary with 
rapidity. 

This requirement might seem counter-intuitive: At the 
critical density separating partonic and hadronic degrees 
of freedom, it is natural to expect (av) to increase as the 
nearly massless color-charged partonic degrees of free- 
dom are liberated. This likely violation, in fact, was the 
original motivation for the derivation of [65| . 

Considering, however, that the weakly interacting sys- 
tem will be far-away from equilibrium until hadroniza- 
tion (we note that the critical system size required for 
statistical hadronization is a topic of considerable con- 
troversy (6|| |67} ) , and that breaking flux tubes emit par- 
tons rather than hadrons throughout the rapidity range of 
the system, it is not unreasonable to suppose (av) main- 
tains its partonic value even in the dilute large rapidity 
regions. 

A more serious problem is modeling the absolute value 
of Vi : as shown in [681 ] , a model of the type described in 
[65| is not able to describe the absolute value of v-i unless 
a was increased to the point where the mean Knudsen 
number is well below unity. In such a regime, the ap- 
proximation of "one collision per particle per lifetime" , 
used to derive Eq. is clearly not appropriate and the 
ideal hydrodynamic limit emerges. 

One can modify the derivation in [65| with the intro- 
duction of a mean field. This way, the magnitude of i>2 
can be acceptably described (provided coalescence holds) 
even with a negligible a (the Vlasov equation limit) (69j . 
provided the lifetime of the system is long enough for 
the mean field description to be appropriate. If this is 
done, and the magnitude of the mean field is proportional 
to dN/dy (a reasonable proposition), one can substitute 
(av) — > (M) in Eq. [T5] ((M), the mean field per parton, 
being an independent constant) and increase vi far above 



the limits set by the K ~ 1 approximation. The scaling 
naturalness inherent in Eq. [T5] is maintained provided 
the system is, at formation, partonic (and described by 
Eq. [12]) throughout all rapidity. 

While such assumptions arc not unreasonable since 
repulsive mean field potentials seem to be suggested 
by lattice QCD [69j], we are far away from assessing 
whether crucial aspects of this scenario (e.g. whether 
(M) oc dN/dy at all energies, and whether it is really 
true that v-i forms at times -C r eq ) fit with QCD [TOl — 
[72j . We therefore limit ourselves with noting that such 
a model is the only one on the market where the strict 
proportionality of V2 with dN/dy, and hence the limiting 
fragmentation of i>2, can be incorporated naturally. 



V. DISCUSSION AND CONCLUSIONS 

The unnaturalness of V2 limiting fragmentation in the 
strongly interacting (hydrodynamic) limit, and the prob- 
lems encountered to model V2 in the weakly interacting 
(Vlasov equation) limit motivate us to look for experi- 
mental ways clarify the situation. 

The LHC might well tell us very quickly if hydrody- 
namic descriptions of vo, are appropriate, since the hydro- 
dynamic prediction [73l - l75l | differs from the extrapolation 
from limiting fragmentation [34], [76| by as much as 40%. 
The preliminary results suggesting that the multiplic- 
ity at mid-rapidity increases faster than dN/dy ~ In \fs 
[77L EH suggest, if the ideas discussed in section III Bl are 
correct, that limiting fragmentation at the LHC will 
break down to the same extent. It will be very inter- 
esting to see if this is the case. 

A promising endeavor for low energy measurements is 
to look for limiting fragmentation in other soft observ- 
ables. In particular, it might be worthwhile to look for a 
similar limiting fragmentation in the average transverse 
momentum (pt) |79j - a break of which could signal the 
"step" |8(| in rapidity space, (pr) ~ T + m n v^, 
and vt should decrease significantly when the initial T 
is around T c due to the dip in c s in the mixed phase 
regime. Similarly, rapidity dependence of R OM t and R S ide 
HBT radii could be used to check whether the scaling 

seen in global R ou t,side ~ \llif) radii [HJ is also local 
in rapidity. Since accounting for R ou t,side a t mid-rapidity 
with hydrodynamics requires adding flow created before 
equilibrium [281 ] , limiting fragmentation of HBT radii will 
pose even greater challenges for hydrodynamics. Such 
measurements will clarify whether the compressibility of 
the system does indeed change with rapidity, as its ex- 
pected to in local thermal equilibrium where it mirrors 
the softness of the equation of state [80j, but not in a 
far-from equilibrium expansion driven by mean fields. 

In this work, we have argued that the experiment ally 
observed scaling of multiplicity with rapidity and t/s [341 ] 
follows from reasonable models of initial partonic dynam- 
ics. Neither the free-streaming limit nor the ideal fluid 
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limit are expected to break it. The situation is however 
different with the equally simple scaling observed with 
V2- It is not clear how this scaling could arise within 
non-ideal hydrodynamics, even if its initial condition mir- 
ror closely the ones that reproduce the scaling observed 
in dN/dy. The goodness of the scaling seen at RHIC, 
therefore, might force us to reconsider weakly interact- 
ing dynamics, where it arises more naturally. Further 
measurements of rapidity dependence of other soft ob- 



servables, within the context of an energy scan, might 
help clarify the situation further. G.T. acknowledges 
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